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O . Abstract 

o 

"^ ' The probabilities for gaps in the eigenvalue spectrum of finite N x N random unitary 

f->i , , ensembles on the unit circle with a singular weight, and the related hermitian ensembles 

on the line with Cauchy weight, are found exactly. The finite cases for exclusion from 
single and double intervals are given in terms of second order second degree ODEs which 
C^ ' are related to certain Painleve-VI transcendents. The scaled cases in the thermodynamic 

limit are again second degree and second order, this time related to Painleve-V transcen- 
dents. Using transformations relating the second degree ODE and transcendent we prove 
an identity for the scaled bulk limit which leads to a simple expression for the spacing 
J^ ' p.d.f. We also relate all the variables appearing in the Fredholm determinant formalism 

JH , to particular Painleve transcendents, in a simple and transparent way, and exhibit their 

scaling behaviour. 



^ 



1 Introduction 

Our intention in this work is to study two intimately related ensembles of random matrices 
with unitary symmetry - the circular Jacobi unitary ensemble consisting of eigenvalues con- 
fined to the unit circle with a spectrum singularity, and the Cauchy ensemble on the real 
line. As we shall see, there is a mapping between these two ensembles which allows us to 
find all the results for the circular Jacobi ensemble from those of the Cauchy ensemble. The 
latter is a classical matrix ensemble on the real line, and as such there is a well established 
method for characterising the gap probabilities in terms of differential equations [|l3|. This 
avoids the clumsy and unnecessary work of recasting all of the formalism onto the unit circle. 
The circular Jacobi ensemble that we consider is one with an algebraic singularity in its spec- 
trum, which in the log-gas picture of the eigenvalue p.d.f. can be interpreted as an impurity 
with variable charge fixed at z = 1 interacting with the mobile unit charges representing the 
eigenvalues. 

In general an ensemble of unitary matrices with unitary symmetry has an eigenvalue 
p.d.f. of the form 

N 

YIW2{ZI) II \Zk-ZjW (1.1) 

1=1 l<j<k<N 

z = e = e^'^*^'^, 9 G [0, 27r), x G [0,L), while an ensemble of Hermitian matrices with 
unitary symmetry has an eigenvalue p.d.f. of the form 

N 

n^2(A,) n \>^k-xj\'. (1.2) 

1=1 l<j<k<N 

The circular Jacobi ensemble with unitary symmetry, denoted as cJUE, is specified by the 
p.d.f. ( |1.1| ) with the weight function 

W2{Z) = \1 - zl'"'' . (1.3) 

In the case a = 1 this is realised by the eigenvalues z = e from the ensemble of (A^ -|- 1) x 
(N + 1) random unitary matrices, with all angles measured from any one eigenvalue. The 
Cauchy unitary ensemble, denoted by CylJE, is specified by ( |1.2| ) with 

^2(A) = (l + A^)-^-^ A G (-CX), oo) . (1.4) 

It results from the cJUE by mapping the eigenvalues from the unit circle to the line tangent 
to its south-most point via a stereographic projection from its north-most point. Specifying 
points on the unit circle via an azimuthal angle 0, so that 6 = corresponds to the north-most 
point and 6 = —tt, vr to the south-most point, this is achieved by the mapping 



A = cot(6'/2) , (1.5) 



and we find 

N 



J]|l-Zi|2« II \zk-Zj\^dei...dBN 

i=l l<j<k<N 



oc 



^ 1 

n (i+A^)^+a n |A.-A,pdAi...dA^. (1.6) 



In (11.21), let us write 



ii;2(A) = e-2^W and 2V'{\) 



g(A) 
/(A) 



(1.7) 



with f{X),g{X) assumed polynomials in A. In O, the defining property of a classical ensemble 
in the random matrix context was identified as the operator 

■f'-g' 



ax 

increasing the degree of the polynomials by one. Thus for a classical ensemble 

degree/ < 2, degree gf < 1 . 

Using this criterium one can check that the weight functions 

e~ , Her mite 

A'^e"'^ (A > 0) , Laguerre 

^(1-A)"(1 + A)^ (-1<A<1), Jacobi 



(1.8) 



W2(A) 



(1.10) 



are all classical. These weight functions all naturally occur as the eigenvalue p.d.f. for certain 
ensembles of Hermitian matrices based on matrices with independent Gaussian elements (see 
e.g. the Introduction section of [P). Furthermore, for each of these cases the probability of 
a single eigenvalue free region which includes an endpoint of the support of the weight has 
been expressed in terms of the solution of a particular non-linear differential equation |13]. 
The same has been done for the Hermite and symmetric Jacobi (a = /3) ensembles when the 
eigenvalue free region consists of two disjoint intervals at either end of the support of the 
spectrum [|l4| . The Cauchy weight function (1^) is significant for the feature of satisfying 
the criterium ( |1.9D and thus being a fourth classical weight function. The orthogonal poly- 
nomial system defined by the Cauchy weight has been investigated in [l^ and 0], where the 



polynomials have been found to be given by certain Jacobi polynomials of pure imaginary ar- 
gument. Consequently we find some similarities with the Jacobi ensemble, but also sufficient 
differences to warrant separate attention. With £'2(0;/; W2(A); A^) denoting the probability 
that there are no eigenvalues in the interval / of an ensemble with eigenvalue p.d.f. ( |1.2| ), the 
probabilities to be calculated are 



S2(0;(s,oo);(1+A2 



-N-a 



■,N), 



(1.11) 



and 



£2(0; (-00, -s) U (s,oo); (l+A^)-'"-''; A^) . 



(1.12) 



The quantity (1.11) gives the probability that there are no eigenvalues in the CyUE greater 
than s, and the quantity in ( |1.12 ) gives the probability that there are no eigenvalues in the 
CyUE with modulus greater than s. In terms of the circular Jacobi ensemble the probabilities 
corresponding to ( 1.11 ) and ( 1.12| ) are 



i?2(0;(0,x);|l- 



-z\'-;N), 



(1.13) 



and 



£2(0;(-x,x);|l- 



|2a. 



N) 



(1.14) 



The details of the Cauchy ensemble are defined in Section 2, including the orthonornial 



polynomials and associated coefficients. In Section 3 the formalism of Tracy and Widom |13] 
giving coupled differential equations for the gap probability and some auxiliary quantities is 
revised. In Section 4 the coupled equations for the single interval are reduced to an ordinary 
differential equation specifying the probability ( |1.11 ), and this equation is shown to be related 
to a Painleve-VI transcendent. This problem and the associated ODE are mapped back to 
the circular Jacobi ensemble and the thermodynamic limit taken of ( 1.13| ), and the resulting 
ODE and its transformation to a Painleve-V type studied. A similar study is undertaken in 



Section 5 for the probabilities of exclusion from the double interval, namely ( 1.12 ) and the 
scaling limit of ( 1.14| ). 



2 The Cauchy Ensemble and Orthogonal Polynomials 

Consider the problems of computing the probabihty that there are no eigenvalues in a region / 
of the spectrum for an ensemble specified by (|1.4|). Suppose I consists of M disjoint intervals, 
so that with the endpoints of these intervals denoted by {aj}^^^, 

M 

1= IJ (02m-l,a2m) • (2.1) 

m>l 

The probability of no eigenvalues being found in this interval is given by the general expression 
(see e.g. §) 

^(0;/) =l + Y^^:^ I dxi... [ dXnPnixi, . . . ,Xn) , (2.2) 

where pn is the n-point distribution function of the eigenvalue p.d.f. For Hermitian matrix 



ensembles with unitary symmetry the eigenvalue p.d.f. is proportional to (1.2), and the 
corresponding n-point distribution function is given in terms of the orthonormal polynomials 
{pj(rc)}j=o,i,2,... associated with the weight function W2{x) according to the formula 

Pn{xi, ... ,Xn)= det [KN{Xi,Xj)]^^.^.^^ , (2.3) 

where 

N-l 

KN{x,y) = [w2{x)w2iy)]^^'^ ^ Pi{x)pi{y) . (2.4) 

1=0 

Substituting into ( p.2| ) , the Predholm theory of integral operators then gives 

E{{);I) =det(I-KAr) , (2.5) 

where K^r is the integral operator with kernel K]\j{x,y) defined on the interval /. A crucial 



point is that (2^) can be summed according to the Christoffel-Darboux formula and so 



written in the special form 



KN[x,y) = , (2.6) 

x-y 



where with ajy denoting the coefficient of x^ in p^ {x 

1/2 



(f){x) = ( -W2{x) ] pn{x) , 

aN J 

(2.7) 

V'(x) = ( -^^W2{x) ] PN-l{x) ■ 
aN J 



Now, the Cauchy weight is defined by ( |1.4D where A'^ G 1j and a G R^. This weight 
defines an orthogonal polynomial system only up to the degree n = N and not for higher 
degrees, however this is of no consequence as we only require those below N. The Cauchy 
weight is closely related to the symmetric Jacobi weight (1 — A^)°^, as can be seen from the 
following integration identity M 



Lemma 1 For analytic functions f{x), such that the following integrals exist and ?R.{a) > — 1, 



(ix/(ix)(l+x^)" = tan(7ra) / da;/(x)(l-x^) 



„2\a 



Using (^), the fact that 



dx(l-x2)°P^"'°)(2;)P^°'°)(x) 



22"+ir2(a+n+l) 



-1 



,(",/3) 



(2a+2n+l)r(n+l)r(2a+n+l) 



(2.9) 



with Pn ' (x) denoting the usual Jacobi polynomial, allows us to conclude that 



p^y{x) = ^-2^+-^ 



n!(iV+a-n-y2)r2(iV+a-n) 
27rr(2iV+2a-n) 



1/2 



P^^"'^-^-'^)(^x) , 



(2.10) 



(",/3), 



are the orthonormal polynomials for the Cauchy weight. From the definition of the Pn ' (x) 



it follows from (|2.10D that the coefficient of x" required in ( |2.7D is such that 



Cy 

Cy 
an 



V2 



n{2N+2a-n) 



{N+a-n+y2){N+a-n-y2) 



1/2 



(2.11) 



A common feature of the classical orthogonal polynomials is that (/>(x) and ^l^{x) satisfy the 
recurrence-differential relations 



m{x)(j)'{x) = A{x)(j){x) + B{x)ip{x) , 
m{x)ip'{x) = —C{x)(j){x) — A{x)il^{x) , 



(2.12) 



where the coefficient functions m{x), A{x), B{x),C{x) are polynomials in x. Making use of 
the differentiation formula 



d 



(2n+a+/3)(l-x2)— P^"'^)(x) = n[a-p-{2n+a+p)x]P^ 



",/9)| 



dx 



X) 



+ 2(n+Q)(n+/3)P^^f (x) , (2.13) 



(",/?) / 



and the three term recurrence for the Jacobi polynomials gives that ( 2.12[) holds for the 
Cauchy ensemble with Q 

m(x) = 1+x , 

^(x) = —ax , 



B{x) 
C{x) 



N{N+2a) 
N{N+2a) 



2a-\ 


2a+l 


2a+l" 


2a-l_ 



1/2 



1/2 



/3o 



70 



(2.14) 



with n = N . 

Tracy and Widom [|l3| have developed a formalism based on the integral operator deter- 



minant formula (^.5|), and which makes essential use of the special structures (2^) and ( 2.12 ) 
to derive a set of coupled differential equations for £/2(0;I) and some auxiliary quantities 
The salient features of this theory will be revised in the next section. 
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3 The General Formalism 



Our goal is to characterise the probabihties ( |1.11| ) and ( 1.12 ) as the solution of certain non- 



linear differential equations. Following ||13|] this is achieved by specifying partial differential 
equations for the quantities qj,pj,u,v,w and R{aj,ak) defined henceforth. 

Definition 1 Let A = A{x,y) denote that the integral operator A has kernel A{x,y). Then 
the kernels p{x, y) and R{x, y) are specified by 

{l-K)-'=p{x,y) , 

(3.1) 
K{l-Ky' = R{x,y) , 

The operator K{1 — K)~^ is called the resolvent and R{x,y) the resolvent kernel. 
Definition 2 For k G Z>o the functions Qt and Pt are defined by 



Qk{x) = dy p{x,y)y^(t){y) , 



Pk{x) = dy p{x,y)y tl){y) , 



(3.2) 



(3.3) 



and their values at the endpoints Oj of I are denoted qkj,Pkj so that 

Qkj = Qk{aj) = lim Qk{x) , 

Pkj = Pk{aj) = lim Pk{x) . 

x^aj 

Where there is no confusion, we denote qoj,Poj by qj,Pj- 
Definition 3 The inner products u, v, w are defined by 

u = {4>\Q) = J dy QoiyMy) , 

V = {^l^\Q) = Jdy QoiyMy) = {<t>\P) = j dy Po{y)Hy) , (3.4) 

w = {ij\P) = j dy PQ{y)i;{y) . 

The coupled differential equations come in two types. There is a set of universal equations 
which are independent of the recurrence-differential equations ( |2.12|) , and a set of equations 
which depend on the details of ( |2.12| ). Let us first present the former. 

Proposition 1 For general functions (j){x),'ip{x) we have the relations 

^logdet(l -/C) = i-iy~'R{aj,aj) , (3.5) 



and for j ^ k, 



R{aj,ak) = , (3.6) 

dj — a-k 



and 



— i?(a„a,; 



'lfR{aj,ak)R{ak,aj) , 



(3.7) 



along with 



for j ^ k anc 



-^ = {-lfR{aj,ak)qk , 
-^ = {-lfR{aj,ak)pk , 



!!-<-'« 



(3,! 



dv 
dak 



{-ifqkPk , 



(3.9) 



dak -^ '^ ^'^ 



The second set of equations, which depend on the details of ( |2.12|) , give the j = k cases 
of ( |3^ and ( |3^ . Now for the Cauchy weight (|L^ we know from ( |2.14| ) the equations ( p. 121) 
hold for 'm(x) a quadratic and A{x),B{x), C{x) linear functions, and thus of the general form 



(3.10) 



m{x) = fj,o + nix + fi2X , 
A{x) = ao + aix , 
B{x) =Po + Pix , 
C{x) = 70 + 71 2^ • 
One then has the following equations [|13| . 

Proposition 2 /n t/ie case that 0(x),'0(x) satisfy the equations ( 2.1^) with coefficient func- 
tions (3.H) we have 

rrii—^ = [ao + aiCj + 7iii - (3iw - H2v]qi 

Otti 

+ [Po + Piai + 2aiu + 2/3if + fi2u]Pi 

2M 



y^^{-lfR{ai,ak)qkmk 



dpi 



(3.11) 



mj— = [-70 - jiOi + 27it; + 2aiw - ii2w\qi 

+ [-ao - aittj + /Jiw - 71 li + ;U2w]pj 



2Af 



^(-l)''i2(ai,aA:)pA;m'fc , 



and 

miR{ai, tti) = [70 + 7iai - 271^ - 2aiw + fj-2w]qf 

+ [Po + PiUi + 2aiu + 2l3iv + H2u]pi 

+ [ao + aitti + 7iu - Piw - ^2v]2qiPi 

+ 2^(-l) nik -— , (3.12) 



and furthermore, 

d 



-— ymiR{ai, ai)] = 2aiqiPi + /3ipf + jiqf 
oai 



where rrii = m{ai). 



2M 



Y,i-'^)''rnkR\a,,ak) , (3.13) 



4 The Single Interval 

We first consider the probability for tlie interval (s, oo) to be free of eigenvalues so that 
ai = s and 02 = cxd. We shall adopt the conventions qi,pi = q,p and R = R{s,s), noting 
that q2,P2 = 0. 

Proposition 3 The coupled differential equations for the finite N CyUE on the interval 
{s, 00) for general a, 13 are 



[hiE2 



w' 
[l + s^)q 
[l + s^)p 



R, 



qp 

2 



-p 



-[as + v]q + [Po - u{2a-l)]p , 
-[70 + w{2a+l)]q + [as + v]p , 
{l + s'^)R = [70 + w{2a+l)]q^ + [I3q - u{2a-l)]p'^ - [as + v]2qp 
[{l + s^)R]' = -2aqp , 



4.1 
4.2 
4.3 
4.4 
4.5 
4.6 
4.7 
4.8 



Proof (sketch) - The first equation follows from ( p.5| ) and (2.5). The next three ( [4.2| 
follow from ( |3.9D and the pair ( 4l5|), ( [4.6D follow from ( p. 11 ). Equation ( [4.7| ) follows from 
(|1|), while ^J) follows from ( pl3|) - □ 

The boundary conditions satisfied by R{s,s) as s — > 00, deduced from the fact that in 
this limit R{s, s) ~ Kj\[{s, s), can be expressed as 



R{s,s) 



DN-I "^ ^ iV!l [a) 2\-N-a-l 



Att T{N+2a) 



X <^ -N{N+2a) 



P: 



{-N-a,-N~a) 



N 



(Z.) 



2 , p(-Af-a -Af-a) / . s p(-Af-a -Af-a) 



2m^sP; 



AT 



z.)p^-r 



'(i.) 



+a^ 



P, 



(-Ar-a,-Ar-a) 



Af-1 



is) 



or considering only the first two leading order terms 

22«r(l+A^+2a)r2(a+l) 



{l + s^)R{s,s) 



2a 



TT r(iV)r2(2a+2) 



2a+l - s~ 



2a+3 



{2N'^+4aN+ia'^+ia+l) 



(4.9) 



(4.10) 



We now indicate how to reduce such a system to a single second order differential equa- 
tion for R = R{s) = R{s, s). 
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Proposition 4 The coupled set of ODEs given in Proposition 3 reduce to the second order 
ODE for a{s) = {l + s^)R{s), 



{l+s^)\a"f + A{l+s^){a'f - Ssa[a'f + Aa\a' - a^) 



+ Sa^saa' + 4 [N{N+2a) - a^s^] {a'f = . 

Proof - We note from ( ^.31 ) and ( |4.8| ) the integral 

{l+s^)R = 2av , 



(4.11) 



(4.12) 



where the constant of integration has been set to zero because v, {l + s'^)R ^ as s — > cxd. 
Now we show how a second integral of the motion can be found. Utilising this first integral 
and (|4.7|) we have 



[70 + w{2a + l)]q'^ + [/3o - u{2a-l)]p'^ - [as + v]2qp - 2av = , (4.13) 

while adding p times ( [4. 51) and q times (^^) yields another equation quadratic in the p and q 

(4.14) 



{l + s^){pqy + [70 + w{2a+l)]q^ - [Pq - u{2a-l)]p^ = . 



Now the idea is to combine these two equations in such a way so that use of the relations ( |4.2| - 



4.4) to express products of q,p in terms of derivatives of u,v,w leads to an exact derivative. 
Such a combination is ( [4.13| ) minus 2a times ( 4.14 ), which when equations (4.2-47^ are 
employed yields 

2a{l + s^)v" + 2asv' + 2vv' - 2av 

+ (2a-l)n'[7o + tt;(2a+l)] - {2a+l)w'[Po - u{2a-l)] = . (4.15) 

This can be rewritten as 

2a{l + s'^)v" + 4asv' 

+ 2vv — 2av — 2asv 

- [Po - n(2a-l)]'[7o + u;(2a+l)] - [70 + w;(2a+l)]'[/3o - ui2a-l)] = , (4.16) 

which is a perfect derivative. Given that the boundary conditions as s — > 00 imply u,v,w ^ 0, 
the integration constant can be determined, leading to the second integral 

[/3o _ n(2a-l)][7o + w;(2a+l)] = /5o7o + {l + s^)a' - sa + ^a^ . (4.17) 

We seek now to utilise all these equations to arrive at a second order ODE for o"(s), and to 
do so we note the pair of equations 



(l + s2)a" = [Po - n(2a-l)]p2 _ [^0 + wi2a + l)]q'' , 
2a 



which is a consequence of (4.14|) and (4.12|), and 



a — 



2a2 



{a + 2a^s)a' = [/?o - u{2a-l)]p'^ + [70 + w{2a + l)]q^ 



(4.18) 



(4.19) 



We seek to eliminate the two terms appearing on the right-hand sides of these last equations 
by forming the difference of their squares and using their cross product from the second 



integral ( 4.17 ). The resulting equation is the ODE stated in the proposition. D 
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4.1 Reduction to Painleve Type 

It is immediately clear that we can make an identification of our ODE with the master 
Painleve equation for second order second degree ordinary differential equations with Painleve 
integrability of Cosgrove and Scoufis, |^] equation (5.1), where in terms of this equation we 
have g{s) = 1 + s^, h{s) = N{N + 2a) — a?s'^, f{s) = 0, and the coefficients ci = 0,C2 = 
1, C3 = 0, C4 = 1, C5 = -a?,CQ = 0, C7 = N{N+2a),cs = 0, cg = 0, cio = 0. 

Proposition 5 The function a{s) = {l + s^)R{s) for the Cauchy Unitary Ensemble on the 
interval (s, oo) is related to a solution of the Painleve VI transcendent uj{t) with the param- 
eters 

a = 1/2 , P = -i/2(iV+a)2 , 7 = i/2(iV+a)2 , <5 = 1/2 - 2a^ , (4.20) 

by the gauge 

t = y^{l-is) , 



V{t) = ^(^(■s)-a^s) , 
and Bdcklund transformations 

' AL0{u}-l){u)-t) \ t(t-l) J ^ Uj{uj-l) '^ \ LO- 

t{t-i) r. uj{oj-i)Y w,. x2 {^-tf 



(4.21) 






(4.22) 



4u;(a;-l) \ t{t-l) J ^^ ^ t(t-l)a;(u;-l) ' 



All that is required to make an explicit relation to one of the Painleve transcendents is to find 
the right gauge transformation to put ( 4.11| ) into the canonical form of (SD-Ia) in [y]. Given 



that our function g{s) has zeros at si = oo,S2 = i,s^ = —i this determines 01,02,03,04 in 
the general gauge transformation 

015 + 02 05(T + 06S + 07 

t = — T) = , (4.23) 

03S + 04 03^ + 04 

leaving 05,00,07. The coefficient 05 is found from matching the cubic terms in the ODE (all 
three cubic terms must be matched with the single coefficient), whilst 05,07 are determined 
by the r]'^,r]rj terms. This determines the gauge transformation as given in ( ^.21 ). The free 
parameters in the canonical form are then found to be 

Ai = 3a2 + iV(iV+2o) , 

(4.24) 
^3 = 3a'^ + 2o2iV(iV+2o) , 

Ai = a'^ + a'^N{N+2a) . 

From these four coefficients one solution for the Painleve- VI transcendent coefficients is given 
in ( 4.2OI ) and the corresponding Backlund transformation by ( 4.22| ). Other solution sets are 



possible, which are manifestations of Schlesinger transformations, and all are tabulated below 
in Table 1. D 
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a = 


/3 = 


7 = 


5 = 


y2 


-y2(A^+a)2 


y2(A^+a)2 


y2 - 2a' 


1/2(1 ±2a)2 


-y2(A^+a)2 


y2(A^+a)2 


y2 


i/2[l ± (iV+2a)]2 


-y2«' 


y2«' 


y2(i-iv2) 


l/2[l±iV]2 


-y2«' 


y2«' 


i/2[i _ {N+2af] 


i/2[l ± {N+a)]^ 





2g2 


y,[l - (N+a)'] 


i/2[l ± (iV+a)]2 


-2a2 





y2[l - (N+a)'] 


i/2[l±a]2 


-y2A^(iV+2a) 


2o2 + i/2Af(7V+2a) 


y2(i-a2) 


i/2[1±g]2 


2a2 - i/2N{N+2a) 


y2iV(iV+2a) 


y2(i-a2) 



Table 1: Painleve parameters for the P-VI transcendents characterising the gap p.d.f. for 
the Cauchy Hermitian random matrix ensemble on the interval (s,oo). 



4.2 Special Cases of low A^ 

In this part we present the calculations for the first two finite- A^ cases, that is A^ = 1 and 
A^ = 2, by direct means using the probability £2(0; I). 



Proposition 6 The probability ( \1-1]\ ) for N = 1 and the associated quantities are given in 
terms of the Gauss hypergeometric function 2-^1 (a, b] c; z) by 



E2{0;I) = X 



a{s) 



r(a+l) {1+s') 



2\-a 



(4.25) 



V^r(a+i/2) E2iO;I) ' 
and similarly for the probability with N = 2 , 

r2(a+l) 



i^2(0; i)=x'- rri^':lM i+s')'^-^x 

Vvrr(a+3/2) 

r(a+2) (l + s^)-"-i 
0Fr(a+3/2) ^2(0;/) 

r(a+i) ^_2,_i 



27rr(a+i/2)r(a+3/2) 



(1+^^) 



2\-2a-l 



a[s) 



l + (2a+l)s2]X + 



r(a+i) 



where 



X 



20Fr(a+3/2; 



V^r(a+i/2) 

s - -2i^i(a + l,a + y2;a + 3/2;-s-2 



sil + s 



2\~a 



(4.26) 



r(a+l) , , 

^2 + ^w J,. s 2F1 a+1, 1/2; 3/2; -s^ 
y/TTT{a+y2) 



(4.27) 



Proof - These follow from the integral representations and transformation formulae of the 
Gauss hypergeometric functions. D 

One can show that these two specific cases are solutions of the second order differential 
equation ( 4.11| ) and satisfy the boundary conditions ( [4.1C ), after noting the differentiation 
formula for the Gauss hypergeometric function. 
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4.3 The Thermodynamic Limit for the Circular Jacobi Ensemble 

Using the stereographic projection we map the extended real hne back onto the circle of 
circumference L via 



vrx 



s = cot -— X G [0, L) 



Now as 



i?2(0;(5,oo);(l+s2)-^-";iV)=exp(-^"dt^) , 
this gap probability is then expressed as 



E2i0;{0,x);\l-e 



2-iTiy/ L 1 2a. 



; N) = exp 



dy a{cot-Ky/L) 



(4.28) 



(4.29) 



(4.30) 



One is then interested in taking the scaling limit, interpreted as the thermodynamic limit, 
in which both N,L ^ oo with a fixed density p = N/L. Evidently the natural quantity to 
consider in this limit is 



— cj(cot TTx/L) \-^ a^npx) 



(4.31) 



and so we investigate the above limit of ( 4.11| ) in terms of this new variable 



Proposition 7 In the thermodynamic limit the scaled variable Ta{x) = xa{x) satisfies the 
second order second degree ordinary differential equation 



2-2 



X (fa) - 4(xfa - Ta)fa - 4a f^ + A{xfa - Ta 



0, 



subject to the boundary conditions 



(y^^)2a+l 



x^o r(a+i/2)r(a+3/2) 



(4.32) 



(4.33) 



Proof - This follows directly from the definition ( [4.31|) under the map ( 4.28| ) of ([4.11J ). D 



Proposition 8 The scaled variable Ta{x) is determined by a solution of the Painleve-V tran- 
scendent u{x) with parameters 



a 



Ml T 2af , (3 = 0, 7 = 0, 5 = 2, 



(4.34) 



and is related by 



1 


xu' 




4u 


u-1 

X 


u' 


4i 


i{u-l) 



a u + 



2 
X U 

(m-1)2 ' 



(4.35) 



n' - (1 =F 2a) 



uiu — l) 



X 



XU 

u-l 
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Proof - The equation ( 4.32 ) falls into the class SD-I.b, equation (5.5) of [^ under r^ i— > — r^ 
and so it is a matter of finding the appropriate gauge transformation. Clearly the parameters 
02 = ^3 = 0, ai = a4 = 05 = 1, and this leaves 06,07. Under such a mapping the canonical 
equation can be compared with our ( [4.32 ), and thus we find ag = 0, 07 = —a^ and Ai = 
—4, A2 = —So?, A^ = 0, A4 = — 4a^. Solving for the Painleve-V transcendent parameters and 
the transformation then leads to the results in the above proposition. Another valid solution 
has the parameters a = 1/2, (i = — 2a^, 7 = 0, 5 = 2. D 

In the case a = the eigenvalue p.d.f. defining the cJUE is identical to the eigenvalue 
p.d.f. for the CUE of random unitary matrices. It is well known that in the scaled thermo- 
dynamic limit the n— point distribution function then has the structure ( p. 3]) with Kj\f{x,y) 
replaced by i^(a;, y) , which in turn has the explicit form 



K{x,y) 



sin TTp{x — y) 



tt{x - y) 
For this scaled ensemble it is a celebrated result of Jimbo et al |jlO| that 



E2(0;(0,x)) =exp 







y 



-dy 



(4.36) 



(4.37) 



where t{x) satisfies the equation ( 4.32 ) with a = 0. 

We have already remarked in the Introduction that the case a = 1 of the cJUE is also 
related to the CUE, by measuring angles from a particular eigenvalue, taken as the origin, 
in the latter. Thus 



^2(0; (0, 9); W2{9) - 1; A^ + l)|given an eigenvalue at 9=0 

But a simple argument shows it is also true that 

N 



^2(0;(0,e);|l-z|2;iV) . (4.38) 



27r 



E2(O;(O,e);«;2(0) = l;iV+l) 



dven an eigenvalue at 0=0 



d_ 

"de 



E2{0;{0,ey,W2{e) = l;N+l) , 



where N/2Tr is the density of eigenvalues at 

2tt d 



N de 



-- 0, so we have the identity 
E2{0;{0,ey,W2{e) = l;N+l)=E2{0;{0,e);\l-z\^;N) 



Putting 



this implies 



^2(0;(0,x);a) 



lim ^2(0;(0,27rx/L);|l-e 



ie\2a. 



■,N) 



(4.39) 



(4.40) 



(4.41) 



-— i?2(0;(0,x);0)=i?2(0;(0,x);l)=exp - / ^^dy], (4.42) 

pdx \ Jo y J 



where the second equality follows from ( [4.30 ) and Proposition 7. 

We have noted the identity (4.42) in a recent Letter |M| presenting the exact Wigner 
surmise form for the spacing probability in the bulk of matrix ensembles with orthogonal or 
unitary symmetries. This is possible (in the case of unitary symmetry) because the spacing 
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probability between consecutive eigenvalues, denoted P2{x) say, is related to -£^2(0; (0, x);0) 
by 



P2{x) = ^E2{0;{0,xy,0) 



where the mean spacing in normalised so that p = 1 • Use of the formula ( [4.42| ) gives 

' n(y) 



Tiinx) 

P2{X) = 2 "3^? 



-dy 



(4.43) 



(4.44) 



X \ Jo y 

where ti satisfies the non-linear equation ( f4.32| ) with a = 1 subject to the boundary condition 

^3 



^^(^^::o3^ 



(4.45) 



Substituting for -^2(0; (0, x); 0) on the left hand side of ( [4.42| ) according to ( ^.SOj ) and the 
result of Proposition 8 shows that ( 4.42| ) implies an identity between transcendents, which 
in fact can be established directly from the theory in Q. 



Proposition 9 The solutions ro(x) and ti(x) specified in Proposition 7 are related by 



Ti{x) = 1 + ro(x) - X 



7"o(a;) 



(4.46) 



Proof - In Proposition 8 for the general Backlund transformations, we have in the case a = 

2 



-To 



'0 



1 

4mo 


xu'q 
uo-1 


no 


X 


11' 


4ii( 


)("o-l) 


"0 



+ 



x'^Uq 



^0(^0-1) 



XUq 
Uq-I 



(4.47) 



so that if we form the quantity on the right-hand side of (4.4(:) then we find 

n2 



1 +To{x) - X 



r'o{x) 

To{x) 



1 

4^0 



xu. 







Mq-I 



Mo 



+ Uq 



x'^Uq 

(«o-l) 



2 ■ 



(4.48) 



But now we note that the solution of the Painleve transcendent for o = 0, uq, corresponds 
to a = 1/2, but is also identical to the one for a = 1, ui, because a is the same (taking the 
negative sign in ( 4.34| )) and thus one can make the identification 



"^^^"^ = 4^ 



XUn 



Uq-I 



Uo 



uo + 



x'^Uq 



(4.49) 



so that ([4.46| ) is established. D 

To understand that ( [4.44 ) is the simplest possible way of representing the spacing prob- 
ability one only needs to observe that the identification of the solutions of two distinct second 
order second degree ODEs with one transcendent can only arise for a = 0, 1 in a = 1/2(1— 20)^^. 
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5 The Double Interval 

We will now consider our second case, which is the probability of eigenvalues being excluded 
from the interval (—00, —s) U (s, 00), in a parallel manner to that of the previous case. We 
make the new conventions oi = —00,02 = —s,a3 = 8,04 = 00, and 93,^3 = q,p, with the 
symmetries q2 = (— 1) y,P2 = (—1) p. We note the parity relation implies v = 0. 

Proposition 10 The coupled differential equations for the finite N CyUE on the interval 

(— 00, — s) U (s, cx)) are 



[lnE2]' = 2R, 

u' = -2g2 , 
w' = -2p^ , 

(l + s2)g' = -asq + [/3o - u{2a-l)]p + ^ ^-q^p 

2('l + s2\ 

(l + s^y = +asp - [70 + w{2a + l)]q - — -qp"^ 



{l + s^)R = [70 + w{2a+l)]q' + [po - u{2a-l)y - 2asqp + 2s(l + s^)i?^ , 

s'~ 



[(l + s')R\' = -2aqp-^^^^qV 



5.1 
5.2 

5.3 

5.4 

5.5 

5.6 

5.7 

5.8 



Proof - These follow in an entirely parallel manner as for the single interval case. D 

Again such a system can be reduced to single second order differential equations for R{s 
or Ro{s). 

Proposition 11 The coupled set of ODEs given in Proposition 10 are equivalent to the 
following second order ODE for a{s) = (l+s^)i?(s), 

S{l-s^)F - 2as + fiilf!! [ah - 2sa' - il+s^)a"] \ 

2{l + s^)a + 2as^{F - as) - s{F - asf 

-is'^iF - asf {N{N+2a) - 2sa - 2as{F - as)} , (5.9) 



where F = ^J a^s^ — 2(1 + s^)(t', with the boundary conditions 

1 r(N+2a+l) _n„ 

s^oo 2^'^+^r(a+i4)r(a+3/2) r(^) 

Alternatively the coupled ODEs can be reduced to the following second order ODE for Rq{s), 
{s{1+s^)^Rq + 2{1+s^){1+2s^)R'q + T [2N{N+2a) + l + s^] - 6aT^ - AT^Y 

= [as - 2{l-s^)Ro]'^ {{l+s^)^{T')^ - 4T^ [T(r+2a) - N{N+2a)]} . (5.11) 
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where T{s) = 2sRo = 2qp, subject to the boundary conditions 

Proof - In order to find our second order ODE we require an integral of the motion (not two 
as for the single interval case because v = 0), and the construction of this parallels the earlier 
case. Firstly we note 

(l + s')M' + [70 + u;(2a+l)]g2 _ [/Jg _ n(2a-l)]/ = , (5.13) 



using ( p7qj5.6| ) and 

0- - [70 + w{2a + l)]q'^ - [/3o - u{2a-l)]p'^ + 2asqp - 2s-^(l + s2)gV = , (5.14) 

which is (|5.7l) rewritten. Now we employ a certain combination of these last two equations, 
namely (|5.14|) + 2a times ( ^.131) , which upon the use of (5^,x4) becomes 



y2{[7o + w^(2a+l)][/?o -n(2a-l)]}' 
+ 2a{l + s^){qpy + Aasqp 

+ a- 2asqp - 2s-^(l + s2)gV = . (5.15) 

However the last two terms of this equation are just sa' and thus the whole left-hand side is 
a perfect derivative. Utilising the boundary conditions as s ^ 00 we arrive at the integral of 
motion, 

[70 + w{2a+l)][Po - n(2a-l)] = /?o7o - 2sa - 4:a{l + s^)qp . (5.16) 

As our starting point in the reduction to an ODE in terms of a we consider the two equations, 
( 15.13 , 5^1^ ), rewritten in the following way 

{l + s^){qpy = [/?o - ui2a-l)]p^ - [70 + w{2a+l)]q^ , 

a + 2asqp-2s-\l + s^)q^p^ = [Po-u{2a-l)]p^ + [jo + w{2a+l)]q^ . (5.17) 

Now we view ( |5.8D as a quadratic equation for qp, with 

.P=^(^[F--], (5.18) 

and use this relation and the definition of F to express the left-hand sides of ( |5.17] ) solely in 
terms of F and a. Finally the difference of the squares of each are taken, and the integral of 



the motion (5.16) is used to eliminate the cross term. One is then left with the first of the 



second order ODEs (|5.9D 



For the second ODE we seek to first eliminate [70 + 'w{2a+l)]q^ and [I3q — u{2a—l)]p'^ by 
employing (^]^) (the sum of these two terms), (p.l3|) (the difference of the two) and (p. 16) 



(the cross product of the two). In this way we find the following quadratic equation for a 
4:ia + T[as-{l-s'^)Ro]\ + 4:T^lT{T+2a) - N{N+2a)\ - {l + s^fT'^ = , (5.19) 



and can then use this equation in conjunction with ( ^.8] ) to eliminate a entirely. The final 
result is the second ODE (|5ll] ). D 



Remark - Knowledge of the above integral ( ^.16 ) allows us to express the auxiliary 



quantities q,p,u,w in terms of the primary ones a,Ro, and this will facilitate further un- 
derstanding of the relationship of these to the transcendents and their scaling behaviour. 
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Proposition 12 The auxiliary quantities q,p are determined by quadrature from the two 
primary quantities a, Rq by the differential equations 



q 2sRq 



2sRo 



+ as-{l+s'^)Ro 



2 ^1/2 

+ 2sa + Aas{l+s'^)Ro - Polo } 



(5.20) 



(l + sY- 
P 



a 



2sRo 



[l+s')Ro 



+ as-{l+s'^)Ro 



2sRo 
subject to the boundary conditions 

1 /2a-l^'/' 



2 ^1/2 

+ 2sa + ias{l+s'^)Ro - /?o7o > 



lis) 



p{s) 



-oo r(a+y2) \2a+l 



1 



2a+l 



1/4 



^- -r(iV + 2a) 1^/^ , , „ 



(5.21) 



s^ocT{a+y2) \2a-lJ { T{N) J 

Proof - One starts by recasting (|5.5D for u in terms of q and the primary quantities, namely 



[/3o-n(2a-l)] 



sRq 



[{1 + s'^)q' + asq - 2{l + s^)qRo\ 



(5.22) 



and similarly for [70 + ti;(2a+l)] by using this result and the integral ( 5.16| ). One can then 
eliminate u, w using these expressions in (|5.17] ) thus relating q, q' solely in terms of a, Rq by 
a quadratic equation, and arriving at ( 5.20| ). The case for p,p' is virtually the same. The 
boundary conditions for q,p follow from the asymptotic forms of (j){s),^lJ{s) as s — > 00 and 
their connection with Jacobi polynomials, (P??, 2.10 ). The sign of the radical is determined 
by a comparison of the boundary conditions for q,p and a, Rq. D. 

The other auxiliary quantities u,w can then be found from quadratures also, given q,p 
using ( ^JsT^ ). 



5.1 Reduction to Painleve Transcendents 

The second degree second order equations of Proposition 11 can also be related to Painleve 
transcendents. 

Proposition 13 The variables Ro{s), F{s), and a(s) are given in terms of the Painleve-VI 
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transcendent uj{x) by the formulae, 



€i2xix-l)uj' -{uj-l){ei{uj + x)-2ax} 
-MKo = r— ^7 -T , [p. 16) 

lyJx{X — 1)UJ 

ei2x{l - x)lo' + ei{u: - l){uj + x) + 2ax 
-iF = —= , (5.24) 

Zy/XUJ 



{2x{x - l)u' - (a; - l)(a; + x)| 



2 



'&^/x^Jo{uJ — l)(a; — x) 

^{u:-l){N{N+2a)uj + a^x] 
iLoyuj — x) 

where €i := ±1 and s = i\fx, and with the parameters, 

« = !' /3=-l{lT2af, 7 = 0, 6 = '^{l - [N+af). (5.26) 

Proof - Under the transformation s i— > —is and Rq i-^ —2iRQ the equation ( ^.llj) is identical 



in form to equation (5.14) in |14], and consequently the rest of the analysis there carries 



through with an adjustment of the transcendent parameters. D 

5.2 Special Cases for low A^ 

In this part we present the analogous results for the first two finite- A^ cases A^ = 1 and N = 2 
by direct calculation of the probability E2{0;I). 



Proposition 14 The probability ( \1.1^ ) for N = 1 and the associated functions are given by 



0rr(a + y2) 

r(a + l) ^_2a-l Z7 /„ , 1 . , w.„ , ./. „-2n 



cr[s 



(l + s2)-'^ 



s-^^-^ 2Fi{a+l, a + 1/2; a + 3/2; -s' 



(5.27) 



F{s) = as + 



2s2Fi(a+l,i/2;3/2;-s2) 

(1 + ^2)-'^ 



S 2^1(0 + 1,1/2; 3/2; -s2) ' 
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and for the corresponding probability with N = 2 is 

I Vvrr(a + 3/2) J 

, , (l + s2)-«-i r 2(a+l)s2 1 



F(s) =as + 



2s [2Fi(a+l,i/2;3/2;-s2)-(l+s2)-'^-i 2i^i(a+2,i/2;3/2;-s2) 

(l + s2)-a-i r 2(a + l)s2 1 



s i2i^i(a+l,y2;3/2;-s2)-(l+s2)-«-i 2^1(0+2,1/2; 3/2; -s2 

(5.28) 



in terms of the Gauss hypergeometric function 2-^1(0, b\ c; z). 

Proof - We proceed in a parallel manner as in the proof of Proposition 6 and the relations 
given there. D 

Again one can show that these two specific cases are solutions of the second order 
differential equations (5. £, 5.11] ) and satisfy the boundary conditions ( 5.10 , 5.12 ), using the 



contiguous and differentiation formulae for the Gauss hypergeometric function. 

5.3 The Thermodynamic Limit for the Circular Jacobi Ensemble 

In Section 4.3 we saw that the stereographic projection ( [4.28[ ) maps, in the limit L -^ 00, 
the interval (s, 00) in the Cauchy ensemble to the interval (0, x) in the scaled circular Jacobi 
ensemble. Similarly this mapping and limiting procedure maps the region (00, —s) U (s,cxd) 
in the Cauchy ensemble to the interval (— x, x) in the scaled circular Jacobi ensemble. Thus 

^2(0; (-00, -s) U {s, cx)); (l+t^)-^-"; N) = ^2(0; {-x, x)- |l-zp"; N) , (5.29) 

and since with a{t) specified by Proposition 11 

a{t) 



(|>00 
-2 J dt-^^. 

we have 



(5.30) 



^2(0;(-x,a;);|l-e2'^*^/-^|2«;Ar) =exp(-^ / dy a {cot Try/ LU . (5.31) 



Furthermore, with 



E2{0;{-x,x);a)= lim E2{0]{-x,x);\l-e^'''y/^f'';N), (5.32) 

N/L=p 
— :0-(cot TTy/L) H^ o-(7rpy) (5.33) 
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( |5.31| ) gives 



S2(0;(-x,x);a) =exp(-2 / dy a{y)] . (5.34) 



It is of relevance to consider the scaled form of the Cauchy kernel under the stereographic 



projection (4.28). Using the relation between the Jacobi polynomial and the hyper geometric 
function 2-^"!, transformations and confluent forms of the latter, and the relationship of iFi 
to the Bessel function one can show Q 

KN{si,S2)ds2 -^ K{xi,X2)dx2 , (5.35) 

where 

Kix,y) = (vrpx)V^(vrp.) V^ Jg^iM^P^Ua-i/^i-Py) - J.+i/2(vrpy) Ja-i/2(vrpx) ^^^^^^ 

2[x y) 

This particular Bessel kernel was first identified in [^]. It has been used in [^ to provide a 
direct derivation of (5.34|), and furthermore specify cr{y) therein as the solution of a particular 



second order ODE. This latter result can be considered as a limiting case of (|5.9| ). We 
summarise here the known results for this kernel, for use later on, and whilst we use the same 
symbols for primary and auxiliary quantities in this case and for the finite Cauchy problem 
the reader should not confuse them. We intend to display the interconnections between the 
scaled Cauchy quantities and the corresponding Bessel kernel ones in the following subsection. 

Proposition 15 The gap probability and associated quantities for the Bessel kernel on the 
interval {—x, x) with 

(p{x) = {x/2f'''ja+i/2{x) , 
i;ix) = ix/2)'/\j,_,/2ix) , 
are defined by the following differential equations with R = R{x, x) 

u = 2q^ , (5.37) 

w' = 2/ , (5.38) 

xq' = xp + {—a + [u — w])q , (5.39) 

xp' = —xq + {a — [u — w])p , (5.40) 

xR = x{q^ + p^) + 2(-a +[u- w])qp + 2{qpf , (5.41) 

(xR)' = q'^ +p'^ . (5.42) 

along with the boundary conditions 






(5.43) 
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Proof - See reference Q, with the mapping x i— > Tipx. D 

For future reference we note the following integral arising from ( 5.37 , 5.38 , 5.38 , 5^40|) , 

2qp = w -u . (5.44) 



Proposition 16 Consider the scaled function cr{y) in ( 5.35 ), and write cri(r) = —2xcr{x) 
r = 2x. Then we have that cti satisfies the ODE 



l\2 



ra'lY + A[-a'-ai + ra[]{{a[) 



-|2 



a^ + CTi — ro"^ 



0, 



subject to the boundary conditions 

0-1 (r) 



-0 r(a+i/2)r(a+3/2)^4 ! 



(5.45) 



(5.46) 



Proof - The result comes from the leading order terms in the ODE ( ^.9D after the change of 
variables (|4.28| ) and introducing the dependent function (|5.33 ). □ 

This is precisely the equation for ai{2x) = —2xR{x, x) that was derived from the coupled 
differential equations (|5.48| - 5^4^ ) in [||. We now want to indicate the precise relationships 
between the scaled Cauchy ensemble quantities a,RQ,q,p and the Bessel kernel quantities, 
which we shall distinguish by affixing the subscript 00, namely the quantities Roo,Qoo,Poo- 
From the arguments in the proposition above we see that 



a{x) = Rooix) , 



after recognising that 



£^2(0; {-x,x);a) = exp -2 



irpx 



dyRooiy) 



Having made this identification we can combine ( 5.41 , 5.42| ) according to 



x—ai-ai = -4a(gooPo 



'^iqooPc 



(5.47) 
(5.48) 

(5.49) 



Then by performing the scaling limit on ( p.SD and comparing with the above equation we 
find a simple scaling for Rq 



NRo{s) = riirpx) , 
along with the identification in the thermodynamic limit 

r{x) = xq{x)p{x) = xqc,o{x)pooix) ■ 



(5.50) 



(5.51) 



Here we define the scaled auxiliary variables by q{s) 1— > q{TTpx) and p{s) >-^ p{Trpx). Utilising 
these two simple scalings we have the following result. 

Proposition 17 The scaled limits of the auxiliary quantities for the Cauchy ensemble on 
the double interval q{x),p{x) are related to the auxiliary variables for the Bessel kernel 
qoo{x),poo{x) by 



p{x) 



2a+l 



r(a+3/2) \2a-l 



1/4 



{V2) 



a+1 



X" 



, ^ / , Poo(^) 



(5.53) 
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Proof- Utilising the scaling forms ( ^.33 , 5?50|) in the expressions for the logarithmic derivatives 
of q in ( 5.20| ) we find 



ax 



-1 



or 



Po 



Poo 

qoo 



(5.54) 



The correct choice can be made and the integration constant found by imposing both the 
boundary conditions (5.21,5.43). The situation for p is entirely analogous. D 



5.4 Reduction of scaled Limit to Painleve Transcendents 

In this part we show how to reduce the auxiliary variables q,p (we now drop the subscript 
cxd) for the Bessel kernel directly to Painleve transcendents. 

Proposition 18 The variable p{x) for the Bessel function kernel with parameter a on the 
interval {—x,x) is related to the Painleve-V transcendent y{x) by 



P 



x 1 + y 
2 1^ 



whose parameter values are 

a = i(l-2a)2, i3 = -±{l-2af 
In terms of this transcendent q{x),ai(2x) are 

/f^[-xy + i(2a-l)(l-y2)] , 

2(1 + 2/) 



7 = 0, 6 



y 



■4y 

1-y 



+ 1/4(1 + y) 



1/40' 



i + y 
i-y 



Proof - We eliminate all variables other than q and p from equations (13)-(17) in 
( |5.44| ), and employ this with the coupled set ( 5.3£ , 5.40 ). Elimination of q using 



leads to the ODE for p 
2p 



P 



[p 



/\2 



2p^ 



xp — ap 
2p2 — X 



-p + -{2p'^ - x) + a{l-a)- 



2p^ — X x(2p2 — x) X x{2p^ — x) 



Then applying the above transformation reduces this to a Painleve-V. D 



(5.55) 
(5.56) 

(5.57) 

(5.58) 
using 

(5.59) 
(5.60) 
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6 Conclusions 

The circular Jacobi ensemble describes a degeneracy of order a in the spectrum of unitary 
random matrices. We have found the probability that the one-sided interval (0, x) to the 
right of the singularity is eigenvalue free, and similarly for the symmetrical interval (— x, x) 
about the singularity. The former probability has been expressed in terms of a Painleve-VI 
transcendent, which in the scaled limit reduces to a Painleve-V transcendent. These calcu- 
lations were performed by first mapping the circular Jacobi ensemble, via a stereographic 
projection, to the Cauchy ensemble of random Hermitian matrices. This latter ensemble is 
classical, having the same essential properties as the more familiar Hermite, Laguerre and 
Jacobi classical ensembles. As the probability of a gap free interval about an endpoint of 
these latter three ensembles are known in terms of Painleve type transcendents from ear- 
lier studies ||l^, 1^, as is the probability for two intervals at the ends of the support in the 
Hermite and symmetric Jacobi cases [|lj] , by studying the Cauchy ensemble we are complet- 
ing a classification program of determining the Painleve transcendents associated with such 
probabilities in the classical ensembles. 

The Cauchy ensemble is closely related to the symmetric Jacobi ensemble. Previous 
studies of the probabilities of gap free regions in this latter ensemble |jl^, |8|, Q have encoun- 
tered third order differential equations, which through various means have been integrated 
to second order equations. In contrast in this study we are able to find an extra integral of 
the coupled differential equations, unknown from previous studies, which provides a direct 
path to the second order differential equation characterising the gap probability. We consider 
that the method of constructing this new integral will apply more generally than just to the 
ensemble treated here and will have important implications for the application of the Tracy 
and Widom formalism to other problems. 

The case a = 1 of the scaled circular Jacobi ensemble has particular relevance to the 
scaled CUE of random unitary matrices, or equivalently to the scaled CUE of random Her- 
mitian matrices. Thus it corresponds to fixing an eigenvalue at the origin in these latter two 
scaled ensembles. The probability of no eigenvalues in the one-sided interval (0, x) of the 
scaled circular Jacobi ensemble is then proportional to the derivative of the probability that 
there are no eigenvalues in the interval (0, x) of the scaled CUE or CUE. This is turn allows 
us to derive the exact Wigner surmise type form ( 4.44| ) for the spacing probability between 



consecutive eigenvalues in the bulk of matrix ensembles with unitary symmetry. 
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